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Abstract. We report on the marginal detection of cosmic shear on sub-arcmin scales with archive data from the 
STIS camera on board HST. For the measurement 121 galaxy fields with a field of view of 51" x 51" are used 
to obtain an rms cosmic shear of ~ 4% with 1.5a significance. This value is consistent with groundbased results 
obtained by other groups on larger scales, and with theoretical predictions for a standard ACDM cosmology. To 
show the suitability of STIS for weak shear measurements we carefully investigated the stability of the PSF. We 
demonstrate that small temporal changes do not affect the cosmic shear measurement by more than ~ 10%. We 
also discuss the influence of various weighting and selection schemes for the galaxy ellipticities. 

Key words. Cosmology: theory, dark matter, gravitational lensing, large-scale structure of the universe 



1. Introduction 

The tidal gravitational field due to the inhomogeneous dis- 
tribution of matter in the Universe causes the distortion 
of light bundles from distant sources. As a consequence, 
the observed images of these sources are deformed. This 
tidal distortion yields a small but observable imprint on 
the distribution of galaxy ellipticities, an effect called 'cos- 
mic shear'. As pointed out by Blandford ct al. (1991), 
Miralda-Escudcj (1991) and |Kaisei| (1992), the statis- 
tical properties of the shear field reflect the statisti- 
cal properties of the (dark) matter distribution in the 
Universe. Refined theoretical predictions, accounting for 
the non-linear evolution of the large-scale structure (e.g. 
Jain k Seljak] 1 997; [Bernardeau et al.| 1997; [Kaisei] 1998; 
Schneider et al. 1998a) and numerical simulations based 



on N-body results for cosmic structure evolution (e.g., 
van Waerbeke et al. 1999, Jain et al. 2000) yield expected 
shear amplitudes of about one percent on scales of a few 
arcminutes, depending on cosmological parameters and, in 
particular, on the normalization of the dark matter power 



spectrum (for recent reviews on the topic, see Melliei 1999; 
Bartelmann & Schneider 2001). 

Because of the small scale of this effect, investigat- 
ing it observ ationa lly is challenging . Early attempts (e.g. 
Mould ct aT] 1994; |Schneider et alj 1998b) where plagued 
by the small data sets. Furthermore, dedicated image anal- 
ysis software was needed to correct for observational ef- 
fccts like atmosph eric seeing, anisotropic PSFs etc. (e.g, 
Bonnet fc Mcllicr 1995; Kaiser et al. 1995 (hereafter re- 
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ferred to as |KSB95| ); [Luppino fc Kaiscij 1997). 

In March 2000, four groups nearly simultaneously 
announced the detection of cosmic shear on angular 
1' and ~ 10' (Bacon et al. 
2000; 
2000) 



scales between 



Kaiser et al 



Wittman et al 



van Waerbeke et al 



2000; 
2000; 



All four groups presented results 
from groundbased wide-field imaging observations, at the 
level expected from currently popular cosmological mod- 
els; particularly noteworthy is the impressive agreement 
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between the four independent results. Maoli et al. (2001) 
measured cosmic shear on an angular scale of a few 
arcminutes, using a survey carried out with the VLT, 
and again found results in agreement with the others. 
Van Waerbeke et al. (2001) presented the currently most 
accurate cosmic shear measurement, based on about 
6.5 square degrees of CFHT wide- field imaging data. 
With the precision of these measurements, interesting 
cosmological constraints can be obtained from these data. 

To extend these results to smaller angular scales where 
the non-linear effects of structure formation arc more 
pronounced, other observing strategies are preferable. 
The noise of cosmic shear measurements is due mainly 
to two sources: intrinsic ellipticity of the source galax- 
ies, and sampling variance. The relative importance of 
these t wo noise terms depends on the angular scale (e.g., 
Kaiser 1998); for small-scale cosmic shear measurements, 
the intrinsic ellipticity becomes dominant. To reduce this 
noise contribution, one should strive for imaging data with 
a large number density of galaxies, i.e., very deep images. 

On the other hand, fainter galaxies have a smaller an- 
gular extent. The observed images are affected by a PSF 
smearing which must be corrected. The smaller galaxies 
need larger corrections and, accordingly, the correspond- 
ing uncertainty is larger. This means that shear measure- 
ments with groundbased telescopes have a natural limit, 
taking deeper images will not provide a significant increase 
in the number density of galaxies which can be used for a 
shear analysis. To measure accurate shapes of very faint 
and small galaxy images, space-based observations are re- 
quired. 

The Hubble Space Telescope currently carries two op- 
tical cameras, WFPC2 and STIS. Data from both of 
these instruments can in principle be used to study cos- 
mic shear. The advantage of WFPC2 is its larger field- 
of-view, whereas the STIS imager has higher through- 
put and smaller pixel size, which is better adapted to 
the diffraction limit of the HST. In addition, it is known 
from weak lens ing studies of clusters with WFPC2 (e.g., 
Hockstra et al. 1998) that the PSF shows relatively large 
anisotropics towards the edges of the three WFC-chips. 
In contrast to wide-field groundbased images, where the 
PSF anisotropy is measured from the stars on the same 
image for which the shear analysis is performed, such a 
procedure is impossible for HST images: the fields are too 
small to detect a sufficient number of stars to model the 
PSF anisotropy. Since there are als o indications th at the 
PSF in WFPC2 changes with time ( [Hockstra et~al| 1998), 
we dismissed the idea of cosmic shear measurements with 
WFPC2 (but see Rhodes et al. 2001). Instead, analysing 



early STIS images, we concluded that the PSF anisotropy 
was very small, and did not change significantly with time. 

With the initiation of a public parallel program with 
STIS in June 1997, we decided to use these archival data 
for a cosmic shear study. The data set that we have anal- 
ysed has been described in Pirzkal et al. (2001; hereafter 
Paper I ) ; briefly, we have produced a set of 498 coadded 



be used to investigate the PSF, whereas others are domi- 
nated by galaxies and can be used for the shear analysis 
(see Paper 1 and Sect. [2] below for details). 

The rest of the paper is organized as follows: In Sect. || 
we briefly describe our data set, the field selection and 
the data analysis. We devote Sect. |^ to a study of the 
PSF anisotropy, using images containing many stars, and 
demonstrate that the variations of the PSF anisotropy in 
time are indeed very small. A description of the cosmic 
shear analysis is given in Sect. [|, including our procedure 
for correcting for PSF anisotropy and smearing. We derive 
the maximum variations of the resulting mean shear per 
galaxy field when using the PSFs obtained from different 
star fields, and conclude that the STIS PSF anisotropy will 
not affect a cosmic shear measurement above the ~ 10% 
level. Results are presented in Sect. [| We obtain a ~ 1.5a 
detection of the rms shear on the scale of the STIS fields, 

with Jl^f) = 3.871^4%- This value extends previous 
cosmic shear measurements to smaller angular scales. The 
influence of various weighting and selection schemes for 
the galaxy ellipticities is studied; in particular, we show 
that using different PSF anisotropy functions for the el- 
lipticity corrections leaves the result basically unchanged. 
We tested our scheme with simulations and verified that it 
is able to recover a cosmic shear variance to within ~ 15% 
of the input value. We conclude in Sect. || with a discus- 
sion of our results and an outlook on the currently running 
Pure Parallel GO Program with STIS. 

2. Catalogue production and field selection 



A preselection of the 498 associations defined in Paper 1 



was done based on a visual inspection of the fields. We 
rejected fields in which the individual exposures were un- 
suitable for our project; reasons to not use a field were: 
bright stars and/or too many stars with diffraction spikes, 
(almost) empty fields (less than 5 objects per field), lower 
image quality than the rest of the images, or failure to 
co nverge o f the iterative co-addition procedure described 



in Paper I. Some examples of fields which were rejected 
are shown in Fig. [l| 

Furthermore we did not use more than one association 
(data taken in a single visit) at a particular pointing, se- 
lecting the one with the highest S/N. We only coadded 
images taken during the same visit for which the same 
guide stars were used so that we could verify that no ro- 
tation exists between successive images. These fields were 
used, however, to test the consistency of the shape mea- 
surement, see Sect. 4.3. 

The fields were analysed using both SExtractor 
(Bertin & Arnouts 1996) and a modified version of the 
IMCAT package ( |KSB95| ; |Erben et al.| 2001) using the fol- 
lowing procedure: 

1. SExtractor catalogues were produced using the param- 
eter file described in Sect. 5.2 of Paper 1 (see also: 



frames, some of which contain many stars and thus can 



http:/ /www.stecf.org/projects/shcar/). In order to se- 
lect well-measurable objects in the catalogues pro- 
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Fig. 1. Some examples of fields which were rejected after a visual inspection: the fields in the top row are dominated 
by one or a few big objects and are otherwise almost empty. In the star fields in the middle row there are many stars 
with diffraction spikes. The last row shows images with bad noise properties and one empty field. 



duced by SExtractor, we applied a 'flag selection', 
disregarding all objects which were flagged internally 
by SExtractor due to problematic deblending and/or 
thresholding. Furthermore, we did a visual inspection 
of all the fields, masking suspect sources like galax- 
ies with clearly separated sub-components and bright 
stars with diffraction spikes and ghosts. 
IMCAT catalogues were produced simultaneously. 
Since IMCAT was developed specifically for shape 
measurements but not for a reliable detection of 
sources, we performed a 'two-code selection': IMCAT 
and SExtractor catalogues were merged requiring an 
object to be detected uniquely by both codes, which 
means that no object had more than one counterpart 
in any of the two catalogues in a radius of 125 mas (5 
subsampled STIS pixels) for the central coordinates. 
The number of objects in the merged catalogue dimin- 
ishes by no more than 5 objects with respect to the 
SExtractor catalogue for the galaxy fields and includes 
98% of objects for the star fields. The objects which 
were rejected were either flagged by IMCAT (near the 
border, negative flux or half-light radii) or there were 
multiple detections by IMCAT where there was only 



one in SExtractor, which is typical for galaxies with 
resolved star forming regions. 

In the final merged catalogue we used size and shape 
parameters from IMCAT, since it was designed specif- 
ically to measure robust ellipticities for faint galaxy 
images and allows for the correction of measured im- 
age ellipticities for shape distortion introduced by the 
PSF (|Luppino fc Kaiser| 1997; [Hockstra et~al| 1998). 



The position and magnitude were estimated with 
SExtractor where we used MAG_ZP=26.38 which is 
the AB magnitude zeropoint for the STIS CLEAR 
filter mode. The objects retained typically have a 
S/N of more than 5. The signal-to- noise ratio is mea- 
sured with the snratio parameter of IMCAT (see 
2001). 



Erben et al 



After catalogue production, we rejected fields with 
fewer than 10 (100) objects for galaxy (star) fields, which 
leaves us with 121 and 51 fields, respectively. 

To select stars, we used size vs. S/N plots in which 
the stars populate a well-defined strip, as can be seen for 
two examples of star fields in Fig. |^. Since the images 
have very different exposure times, we decided to use S/N 
rather than magnitude to be able to use common criteria 
for all fields. For the rest of the paper we assume that stars 
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Fig. 2. Size vs. S/N diagrams for the star fields 
o3zf01010_3_ass (left) and o48b41010^_ass (right). The 
half-light radius rh is given in subsampled STIS pixels. 
One can clearly distinguish the strip populated by stars 
around rh = 2.3 (57.5mas). The strip at w0.8 pixels seen 
in the left p anel is due to single noisy pixels (see also Fig. 6 
of Paper I). 
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Fig. 3. Mean of the ellipticity components e\ (left) and e 2 
(right) of the star fields vs. exposure start time (Modified 
Julian Date). The straight lines show the mean over all 
the fields, the dashed lines show the la dispersion. The 
circles show the mean over stars in bins between 5.06 x 10 , 
5.08 x 10 4 , 5.10 x 10 4 , 5.12 x 10 4 (MJD) with the error 
bars showing the la dispersion. 



have 2.1 < < 2.6 in subsampled STIS pixels (52mas < 
r h < 65mas) and S/N > 10; galaxies are selected with 
rh > 2.6. The half-light radius is measured with IMCAT 
and is the radius at which half of the flux of an object is 
included. 

The mean number of selected galaxies per association 
is 18 on the galaxy fields (25/arcmin 2 ), which is similar to 



of filter scales and using the scale which gives maximum 



the galaxy number density found by Rhodes et aJJ (2001) 
in the 'Groth Strip' with I < 26mag; however, the way 
galaxies are selected in their paper is different from ours. 

3. Analysis of PSF anisotropy 

Our main scientific motivation for using data from the 
STIS Parallel Survey is their suitability for the detection 
and measurement of cosmic shear. Since the expected dis- 
tortion of image cllipticities on the STIS angular scale is a 
few percent, any instrumental distortion and other causes 
of PSF anisotropy need to be understood and controlled 
to an accuracy of < 1%. In particular, the PSF anisotropy 
needs to be either small in amplitude or stable in time. 
The shape of the PSF anisotropy is estimated using the 



KSB95 ellipticity parameters. These describe the elonga- 
tion and orientation of the equivalent ellipse which best re- 
produces the PSF shape, as determined from the weighted 
second order brightness moments, 



Qtj = / d 2 ew(9)d i e j f(e), 



(i) 



where 9 = {9 1 ,6 2 ),9 = 
tive to the centroid 6 C 



\6\, and angles are measured rela- 
of the object, as defined by 



d 2 ew(e c )ej{e c ) = 0. 



(2) 



f(9) is the observed surface brightness and W(0) is a 
weight function, which we take to be a Gaussian with an 
appropriate filter scale. This scale is obtained by deter- 
mining the significance of detecting a peak for a range 



significance as described in |KSB95 . 

The complex ellipticities are then defined as 

Qn - Q22 + 2iQi2 
g = ■ 

Qn + Q22 

for an object with elliptical isophotes 

1-r 2 



. 2it> 



1 



(3) 



(4) 



where r = b/a is the axis ratio of the corresponding ellipse 
and its position angle is $ = 0.5 arctane2/ei. 

In Fig. I the mean values of the two ellipticity compo- 
nents of stars over the whole field are shown as a function 
of the exposure date. The mean ellipticity of all fields is 
1.5% for e\ and 0.5% for e 2 with a dispersion of 1%, which 
is sufficiently small for our cosmic shear analysis, as will 



be shown in Sects. 4.2 and 5.3. Only a small fraction of the 



scatter is due to the different distribution of stars over the 
fields. If we divide the star fields into time intervals, the 
mean ellipticities agree with each other on the la level, 
therefore the anisotropy can be considered to be constant 
over the time period covered. 

In addition to the variation from field to field (i.e., in 
time), we also find a spatial variation of the PSF within in- 
dividual fields. This effect is shown for two fields in Figs. |J 
and ||, bottom left panels. We fit the ellipticities with a 
second order polynomial at the position 8 = (x,y) on the 
CCD: 

e a {9) = a a0 + a a ix + a a2 x 2 + a a3 y + a ai y 2 + a a5 xy, 

(5) 

where a — 1,2. For the fitting an iterative procedure was 
used: for stars selected by rh and S/N, as described in 
Sect. |2| a first guess of the polynomial was calculated by 
solving the normal equations. From this initial fit, a dis- 
persion was calculated by 

2 Si( e t<(^) ~ e ai) fc\ 
°<* = ~1\T AT ' ^ 
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Fig. 4. For the star field o3zf01010_3_ass the distribution 
of the ellipticities of stars are shown before (top left) and 
after (top right) correcting for PSF anisotropy. The bot- 
tom left panel shows the spatial distribution of the ellip- 
ticities across the STIS field, the bottom right the fitted 
second-order polynomial at the star positions. The length 
of the sticks indicates the modulus of the ellipticity, the 
orientation gives the position angle. The mean ellipticities 
before and after the correction are given at the bottom. 
Note that the mean ellipticity after the correction is zero 
and that the dispersion decreases. 

Only stars for which the ellipticity components did not 
deviate more than 3tr in both e\ and e-i from the first fit 
were used for the second fit. The 3cr clipping was then 
applied to all the preselected stars (not only the left-overs 
from the first fit) to avoid biasing. This procedure was 
repeated until it converged, which was typically after no 
more than 10 iterations, with 5-10% outliers disregarded. 
For the two star fields in Fig. ||[ the ellipticities and the 
values for the fitted polynomials at the star positions are 
shown in Figs. || and [|. 

For the anisotropy correction of the galaxies each 
star field polynomial is applied to each galaxy field [see 
Sect. 4.2]. The time period when the galaxy fields were 
observed is well covered by the observation times of the 
star fields. 

To check the stability of our coaddition procedure, we 
also analysed individual exposures with respect to PSF 
size and anisotropy. The stars in the individual exposures 
were selected by l.K r h < 1.4 STIS pixels and S/N > 10. 
We find that the stars in the drizzled images are slightly 
larger («0.1 subsampled STIS pixels) than in the individ- 
ual exposures [see Table |l|. The broadening of the PSF 
is a known property of drizzle if one uses PIXFRAC=1 




500 1000 1500 500 1000 1500 

X X 

(el) = 0. 0211-0.020 {e2) = 0.000±0.0 14 (e 1')= - 0.000 ±0.0 1 B (e2') = -0.000±0.014 

Fig. 5. Same as Fig. |]for the star field o48b41010_3_ass. 



simulations carried out to check our analysis of the cos- 
mic shear; in the simulated star fields we find the same 
increase in sizes of stars (« 5%) from the individual im- 
ages to the drizzled ones, as we observed in the archive 
data. 

From the results in Table [l] we find more stars in the 
drizzled images than in the undrizzled ones. For associa- 
tions where many individual exposures have been coadded 
this can be naturally attributed to a gain in S/N. However, 
we also find more objects if just one individual image is 
drizzled, as can be seen in the last two rows of Table for 
the star field o4xcll010_l_ass. This is a result of the com- 
bination of drizzling and IMCAT: an object which has a 
size of only one STIS pixel or slightly larger is distributed 
over more subsampled STIS pixels in the drizzling pro- 
cess. Therefore an object which was rejected by IMCAT 
on the undrizzled image because of its size of only w 1 
pixel (and therefore high probability to be noise) is now 
detected significantly. 

For the individual undrizzled images of some of the 
star fields, we also analysed the PSF anisotropy and we 
find very short timescale variations of the anisotropy pat- 
tern over the fields, as can be seen in Fig. ^ for two in- 
dividual exposures of the star field o3zf01010_3_ass. The 
short timescale variations of the PSF anisotropy pat- 
tern are likely to be due to "breathing" of the telescope. 
However, they do not affect our analysis, as will be shown 
in Sect. 



4.2 



(see Fruchter & Hook 1998). In Sect. 5.4 we describe some 



4. Shear analysis 
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Table 1. Mean half-light radius and error on the mean for 
individual exposures (in STIS pixels) and for the corre- 
sponding coadded association (in subsampled STIS pixels) 
for all stars on an image. The first Col. identifies the in- 
dividual exposure or the coadded image. The second Col. 
yields the number of stars used to calculate the mean. The 
third and fourth Cols, show the mean half-light radius and 
the error on the mean for individual exposures in STIS 
pixels (50 mas) and for the coadded images in subsam- 
pled STIS pixels (25 mas). The last two Cols, show the 
date of observation idatc (MJD) and the exposure time 
t cxp in seconds. 
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4.1. Cosmic shear theory 

The rms shear in a circular aperture with angular radius 
is related to the power spectrum of the surface mass 
density k by 



/■oc 

(f)=2ir dssP K (s)[I(s6 
Jo 



(7) 



where I(n) :— Ji(n)/(7r?7) and Ji is the Bessel 
function of the first kind. P K is in turn related 
to the three-dimensional po wer spectrum by a 
simple project ion 
Kaiser| 1992, 



Bartclmann fc Schneider 2 001) 



Blandford et al 



1991, 
1999, 

A similar formula for a 



Schneider et al. 1998a, Mellier 



square aperture is given by Kaiser (1992). 

Eq. (Q) is correct if one measures the shear in a circular 
aperture. To quickly compare our result from the square 
STIS field to theoretical predictions, we calculate the the- 
oretically expected value for a circular field with the same 
area, which would then have an effective radius of about 
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Fig. 6. We show the values of the polynomials for the 
anisotropy correction of the two individual (undrizzled) 
exposures o3zf01040 (left) and o3zf01090 (right) which 
are members of the association o3zf01010_3_ass shown in 
Fig. |]. The two individual fields were taken with a time 
difference of only 30 minutes and demonstrate the very 
short timescale variations of the PSF anisotropy pattern. 

30". The error from this approximation is much smaller 
than the statistical error bars on our present result. 

Let ej n denote the complex ellipticity of the i-th galaxy 
on the n-th field, then the quantity we measure for each 
field is 

1 ^ 



7r, 



N„(N n - 1) ^ 



(8) 



where N n is the number of galaxies in the n-th field, or 



In 



(9) 



where Wi n is the weight of the i-th galaxy in the n-th 
field, and e* denotes the complex conjugate. This is an 
unbiased estimate of the cosmic shear dispersion in the 
n-th field. (Note that 7^ is not positive definite.) From 
this, one obtains an unbiased estimate of the cosmic shear 
dispersion: 



1 



where Nf is the number of galaxy fields, or 



(10) 



(11) 



where we weight each field by the number of galaxies per 
field to minimize Poisson noise. Since our galaxy fields 
have a large spread in exposure times, and therefore in 
number of galaxies per field, we use Eq. (11) for our anal- 
ysis. 

Assuming Nf galaxy fields with the same number 
N n = N of galaxies per field and the same redshift distri- 
bution of the sources, the errors from the intrinsic ellip- 
ticity distribution and the cosmic variance in the absence 
of kurtosis are given by 

r 2 



V2- 



(12) 
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and 



)/^Vf- 



(13) 



With an rms intrinsic ellipticity <r s m \/2 x 26% [see 
Fig. 11 1 and an expected shear of a few percent on the 
STIS angular scale we find of ntI . ps 10~ 3 and ct^ v ps 10~ 4 , 
showing that the noise from the intrinsic ellipticity distri- 
bution dominates over the sampling variance. 

4.2. Anisotropy correction 

As shown in Sect. ||, the STIS PSF is remarkably sym- 
metric, i.e., it exhibits a very small (< 1%) anisotropy 
component. Furthermore, it does not vary significantly in 
time (a « 1%), as shown in Fig. |[ 

Of the 51 star fields, 21 were selected for having a 
good spatial coverage and a small intrinsic dispersion in 
the ellipticities of stars, which allowed us to obtain good 
fits to the anisotropy pattern, to minimize the noise in the 
PSF correction. 

The galaxy ellipticities were corrected for PSF 



anisotropy according to KSB95 as follows 



The total response of a galaxy ellipticity to a shear and 
the PSF is given by 



P 7 7 + P Bm q* 



(14) 



where e and e s are the observed and intrinsic ellipticities, 
respectively, and 



pi = p 



sh 



PS, 



(15) 



where the tensors P s h and P sm can be calculated from the 
galaxy light profile with the same filter scale of the galaxy 
as for the weight function in Eq. ([!]); see KSB95. The ratio 
(Psh/Psm)* = (trP* h /trP* m ) is calculated from stars with 
the filter scale of the galaxy. The second term in Eq. ( |l5| ) 
accounts for the circular smearing of the isotropic part of 
the PSF. The stellar anisotropy kernel q* which is needed 
to correct for PSF anisotropy can be calculated by noting 
that for stars e* = and 7* = 0, so that 



q* = (PL) 



-V. 



(16) 



The anisotropy-corrected ellipticity is then calculated by 



e - Psm?* 



(17) 



In Fig. ^ we show for one galaxy field the uncorrected 
ellipticities of the galaxies in the left panel and the mean 
of the anisotropy-corrected ellipticities over all galaxies in 
the field when correcting with different star field polyno- 
mials in the right panel. There is a shift in the negative 
ei direction compared to the mean uncorrected ellipticity, 
as is expected [see Sect. [|. The dispersion between dif- 
ferent corrections is much less than 1% which shows that 
variations of the PSF anisotropy are small. 

In Fig. U the mean ellipticity for all the galaxy fields is 
shown with and without the PSF anisotropy correction. It 




-0.016 -0.014 -0.012 

<e™> 



Fig. 7. For the galaxy field o46p02010.3.ass: in the left 
panel the uncorrected ellipticities of individual galaxies 
are shown, in the right panel the mean of the anisotropy- 
corrected ellipticity (e am ) over the galaxy field for cor- 
rections with different star fields; the square shows the 
mean over the uncorrected ellipticities. Note the different 
scalings of the panels. 



illustrates that the PSF anisotropy correction changes the 
mean ellipticity by an amount typically smaller than 1%. 
Also, the dispersion between different PSF models from 
different star fields is much less than 1%, which means 
that the changes of the PSF anisotropy seen in different 
star fields are sufficiently small to allow us to use one (or 
a suitable combination) of them for the actual analysis. 

The effect of PSF anisotropy corrections is consider- 
ably less than the expected cosmic shear signal, which 
confirms our expectation that the image quality of STIS 
is very well suited for our project. We also note that the 
mean PSF anisotropy correction points towards the neg- 
ative ei direction, in agreement with the stellar ellipticity 
plotted in Fig. 0. 



4.3. Smearing correction 

The smearing corrected ellipticity of each galaxy is calcu- 
lated by 



(18) 



see eqs. (|14|) and (|17|), which is an unbiased (provided 
that (e s ) = 0) but very noisy estimate of the shear 7. 
We apply a scalar inversion of (P 7 ) -1 = 2/trP 7 since it 
is less noisy than the full tensor inversion, as found by 
|Erben et al] (2001). 

To calculate P 7 [see Eq. (pj|) 1 for each galaxy, we need 
to estimate (Psh/Psm)* from the light profile of the stars 
with the filter scale of the galaxy [see Eq. ([!])] . We find that 
(P s h/Psm)* is spatially constant, as for most gound-based 
telescopes (Erben et al. 2001) and calculate the mean of 
this quantity over the stars in a given field for different fil- 
ter scales. We find that (Psh/Psm)* varies from star field to 
star field (i.e., in time), as does the PSF anisotropy If we 
take the mean over the star fields for different filter scales, 
we find that (Psh/Psm)* increases with filter scale and is 
constant for large objects [see Fig. |). This behaviour is 
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Fig. 8. For 121 galaxy fields we plot the mean uncor- 
rected ellipticity of galaxies (triangles) as well as the mean 
anisotropy corrected ellipticity (squares). The error bars 
attached to the squares denote 3 times the dispersion of 
the field-averaged corrected ellipticities when the differ- 
ent PSF model fits are used; the error on the mean is 
much smaller than the symbols used. The shift of the cor- 
rected mean ellipticities towards negative e\ is expected 
from the behaviour of the stellar ellipticities plotted in 
Fig. ||. The big triangle and big square in the centre de- 
note the mean over all galaxy fields of the uncorrected 
and corrected mean ellipticities, respectively; the size of 
the symbols represent the la errors on the mean. The 
circle shows the origin for reference. The mean shear is 
compatible with zero. 

expected since small objects are typically more affected 
by smearing of the PSF and therefore have a larger P* m . 

The error bars on (P s h/Psm)* for different filter scales 
are due to a variation in the size of the stars from field to 
field, as is shown in Fig. [w| where we plot the mean half- 
light radius of all the stars in a star field vs. the time when 
the field was observed. The mean size of stars varies by 
about 0.2 subsampled STIS pixels. We find that in fields 
with larger mean half-light radius, (P s h/Psm)* is bigger for 
all filter scales, i.e., the smearing correction is larger. The 
increasing error bars for larger filter scale are due to the 
fact that stars are small («2.3 subsampled STIS pixels) 
and therefore we detect more and more (background) noise 
if we go to larger filter scales. 

Since for the smearing correction one has to divide by 
P 7 [see Eq. (p!s|)], objects with small values of P 7 can get 
unphysically big ellipticities, which dominate the cosmic 
shear signal even after introducing a weighting of galaxies 
as shown in the following section. We therefore decide to 
introduce a cut in P 7 in addition to the weighting, requir- 



5 10 15 

filter scale (pixel) 

Fig. 9. The mean of (Psh/Psm)* over all stars which were 
used for the anisotropy correction is shown for different 
filter scales. The error bars show the dispersion between 
different star fields. 



ing that 'good' galaxies should have V det P 7 > 0.2 (in 
short: P 7 > 0.2). The effect of different cuts in P 7 on the 
cosmic shear measurement is discussed in Sect. 5.2. 




5.06x10 4 5.08x10* 5.1x10 4 

ExpStartTime (MJD) 

Fig. 10. The mean half-light radius of all stars in a star 
field vs. the time when the star field was observed. The 
filled squares represent those star fields which were se- 
lected for the anisotropy correction of the galaxy fields. 
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Fig. 11. Histogram of the probability distribution of the 
fully corrected ellipticities of all galaxies with P 1 > 0.2 in 
all fields, ei solid, e 2 dashed. The curves show Gaussians 
with mean and dispersion from the histograms. The mean 
ellipticity is very close to zero and a\= <j% = 26%. 

In Fig. [ll] we show the probability distribution of the 
fully corrected ellipticities for all galaxies with P 7 > 0.2 
on all galaxy fields. The mean ellipticity is compati- 
ble with zero and the dispersion in both components is 
26%. The dis persion is consistent with that found by 
~ (1998), 



Hudson et al 



if one takes into account the dif- 
ferent mean magnitude of the galaxies and the different 
pixel scale. The kurtosis is consistent with a Gaussian; the 
skewness is slightly larger than one would expect given the 
number of galaxies. 

For the associations which were observed several times 
at different visits, we compared the catalogues of ob- 
jects obtained from the different exposures and find that 
the number of reliable detections (S/N > 2) differs. 
Comparing the raw as well as the fully corrected ellip- 
ticities, we find that they agree very well regarding orien- 
tation angle and modulus of the ellipticity, with the ex- 
ception of a few faint extended objects with low signal-to- 
noise ratio (S/N < 10). 

4.4. Weighting scheme 

The PSF correction we apply to the galaxy ellipticities am- 
plifies the errors in the measurement, which can eventually 
produce unphysical ellipticities much larger than unity, as 
noted in the previous section. For weak shear the ensem- 
ble of corrected ellipticities measured is the best estimate 
one has of the intrinsic (unlensed) ellipticity distribution. 
But this distribution is broadened and the high ellipticity 
tail is more pronounced because of the noise amplification. 



f 




f 
fi 


















Fig. 12. In the left panel we show the dispersion ctnn of 
individual galaxy ellipticities for the 20 next neighbours 
in the r^—S/N parameter space vs. half-light radius. In 
the right panel the weight deduced from this dispersion is 
shown. 



With the weighting scheme we want to minimize the influ- 
ence of objects with high ellipticities which most probably 
originate from noise. 

Our expectation is that objects which are small and/or 
have low S/N ar most affected by noise. For each galaxy 
we search for the 20 next neighbours in a two dimensional 
parameter space of ha lf-light radius and signal-to-noise 
ratio (see Erben et al. 2001) containing all galaxies from 
all galaxy fields. Since the average number density on our 
galaxy fields is only 18, we have to put all galaxies together 
in one catalogue to do the next neighbour search. For the 
STIS images the galaxies have very similar properties in 
the chosen parameter space, therefore combining them in 
one catalogue is maintainable. 

Since the scalings in and S/N are very different we 
transform the coordinate axes in the following way: we 
sort the objects in ascending order in both coordinates 
separately and assign the running number as the new co- 
ordinate. The 20 next neighbours are then found in this 
new parameter space. For each galaxy we calculate the 
ellipticity dispersion from the M next neighbours 



'NN 



Al 



1 M 



(19) 



which we assign as the measurement error to the galaxy. 
The weights for the galaxies are obtained by a simple 
wnn = 1/°nn weighting 

The dispersion onn and the corresponding weights 
kjnn are shown in Fig. The correlations seen in the plot 
are due to objects with similar properties, e.g. at sa 7 
there are several elliptical galaxies with a compact, bright 
centre. 



5. Results 

5.1. Estimation of the cosmic shear 

The procedure to estimate the cosmic shear is the follow- 
ing: 
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PSF correction: The raw ellipticities of the galaxies are 
corrected for PSF anisotropy with the polynomials from 
the 21 star fields. For each of the star fields we also calcu- 
late (Psh/Psm)* for different filter scales and from that 
P 7 for the filter scale of the galaxy. For each galaxy 
we then have 21 anisotropy corrected ellipticities and 21 
smearing corrected tensors P 7 . The mean correction is 
obtained by then averaging over the star fields. From the 
mean anisotropy-corrected ellipticities and the mean P 7 
for each galaxy we calculate the fully corrected ellipticity 
[see Eq. @]. 

Cosmic Shear: The fully corrected ellipticity is an un- 
biased (although very noisy) estimate of the shear. 
Therefore, we use it to calculate the cosmic shear esti- 
mator for each galaxy field, given in eqs. (^) and @. By 
taking the average over all the galaxy fields as in Eq. ( pi] ) 
we obtain the measured cosmic shear value in our data. 
The error on the cosmic shear result is calculated as the 
statistical error on the mean, which for a weighted mean 



(20) 



is given by 



E^ 2 



(21) 



We also compute the confidence level of the results 
from the dispersion of the probability distribution of the 
cosmic shear estimator in Eq. (|ll) obtained by randomiz- 
ing the galaxy orientations [Fig. 13[] . If we use all galaxy 
fields and weight individual galaxies according to wnn we 
find in 2.9% of the randomizations a value higher than the 
actually measured one. 

Our results for the cosmic shear and the corresponding 
errors [see Eq. (EH])] are shown in Table || for all galaxy 
fields and for fields with more than 10 (15) galaxies per 
field. Note that in some of the fields there are fewer than 10 
galaxies per field (but still more than 5), since we restrict 
our analysis to "good" galaxies, i.e., objects with > 
2.6 subsampled STIS pixels [see Sect. § and P 7 > 0.2 
[see Sect. L3 . Using only fields with a higher number 



density of galaxies we find a larger cosmic shear signal, 
whereas the error stays more or less the same although we 
average over fewer galaxy fields, even if we do not apply 
any weighting to the galaxy fields. 

We also investigated the behaviour of the cosmic shear 
estimator in bins of roughly equal number of fields for 
different number densities of galaxies and find (j 2 ) — 
(-8.1 ± 26.5) x 10~ 4 for fields with 1< N < 10, (j 2 ) = 
(12.0 ± 18.11) x 10~ 4 for fields with 10 < TV < 16 and 
(7 2 ) = (23.4 ± 14.0) x 10~ 4 for fields with N > 16 when 
weighting individual galaxies with u>nn and galaxy fields 
with W{ — N. This shows clearly that the cosmic shear 
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0.004 



Fig. 13. Probability distribution of the cosmic shear es- 
timator calculated from the data by randomizing the ori- 
entations of galaxies with P 7 > 0.2 in all galaxy fields. 
Individual galaxies are weighted with wnn and the galaxy 
fields are weighted by the number of galaxies on each field. 
The vertical line indicates the actually measured value [see 
Table |; 1st row]. 



estimate increases if we use fields with higher number den- 
sity of galaxies. These fields have typically a larger expo- 
sure time and therefore probably probe higher redshifts. 
The main source for the better signal in these fields is the 
reduced error which is due to a better sampling of the 
intrinsic ellipticity distribution. 

5.2. Effect of selection of galaxies and weighting 

In the first three rows of Table || we present the results 
for the cosmic shear for different cuts in P 7 , weighting 
individual galaxies with u>nn- If we lower the cut in P 7 , 
including galaxies whose high ellipticity is most probably 
only due to noise, we find that the cosmic shear estimator 
and the error increase. 

In the second block of Table || we compare the cosmic 
shear estimator for weighting or not weighting individual 
galaxies. We find that the signal increases by one half for 
all bins if we apply weighting, whereas the error increases 
only slightly. 

In the last block of Table || we give the values for 
our cosmic shear estimator applying different weights to 
the galaxy fields. We compare the results for applying no 
weighting at all, weighting by Wf — N, which is the opti- 
mum weight if the error is only due to Poisson noise and 
weighting by Wf = N 2 . We do not discuss a weighting 
by exposure time, which would give deeper fields with a 
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Table 2. Results for the cosmic shear estimator and errors for different minimum number of galaxies per field. N is 
the number of galaxies per field, Nf is the number of galaxy fields with N > N m i n . The first block shows the results 
for different cuts in P 7 , where we weight individual galaxies with w = wnn and the galaxy fields with Wf = N. In 
the next block the results are given for weighting individual galaxies or not. The last block shows the results for the 
cosmic shear if we apply different weights to galaxy fields, weighting them equally (Wf = 1), weighting by the number 
of galaxies per field (Wf = N) or by the square of the number of galaxies (Wf = TV 2 ). Note that we repeat the result 
with P 7 >0.2, to = wnn, Wf = N in each block for easier comparison of the results. 





all 

Nt (f) a ( _ 2> 
x 10 4 x 10 4 


N > 10 
xlO 4 xlO 4 


N > 15 
N t (7 2 > ^- 2) 
xlO 4 xlO 4 


different cuts in P 7 ; w — wnn; Wf = N 


P 7 >0.2 
P 7 >0.1 
P 7 >0.0 


121 14.96 11.61 
121 20.14 12.37 
121 20.41 12.34 


90 16.90 11.42 
92 21.39 12.02 
92 21.49 12.08 


49 26.33 12.28 

50 28.49 12.46 
54 31.96 14.05 


different weighting of individual galaxies; P 7 >0.2; Wi = N 


w =1 

W = 1UNN 


121 9.18 10.17 
121 14.96 11.61 


90 10.29 10.05 
90 16.90 11.42 


49 19.03 10.68 
49 26.33 12.28 


different weighting of galaxy fields; P 7 >0.2; w = ionn 


W f = 1 
Wt = N 
W f = N 2 


121 10.68 11.81 
121 14.96 11.61 
121 15.32 16.74 


90 14.73 11.26 
90 16.90 11.42 
90 15.77 16.27 


49 31.29 12.36 
49 26.33 12.28 
49 19.61 16.08 



higher cosmic shear signal a larger weight, since our er- 
rors are mainly due to Poisson noise. Of course there is 
a strong correlation between exposure time and number 
density and therefore the weighting by N reflects the deep- 
ness of the fields, but there is also the effect of clustering. 
If we increase the weight of fields with many galaxies per 
field, the cosmic shear estimator increases if we take all 
fields, increases and then decreases if we take only fields 
with more than 10 galaxies and decreases if we take more 
than 15 galaxies per field. The error on the cosmic shear 
is the same for the Wf = 1 and Wf = N weighting and is 
much larger for the Wf = iV 2 weighting. This weighting 
obviously gives a too large weight to fields with a very 
high number density of galaxies. 

5.3. Effect of PSF corrections 

In order to estimate how much the PSF corrections affect 
our cosmic shear result we did the cosmic shear analysis 
without either anisotropy or smearing correction or both. 
The corresponding results are found in the first block of 
Table ||. Without the anisotropy correction but includ- 
ing the smearing correction (e raw ,P 7 ) the result does not 
change very much, which confirms that the STIS PSF 
anisotropy is sufficiently small. However, if we leave out 
the smearing correction (e am ,P s h and e raw ,P s h), we find 
that the cosmic shear estimator decreases by a factor of 
roughly three and also the dispersion is much smaller. 
This is not an effect of the different number of galaxy 
fields included in the averaging. The cosmic shear esti- 
mate is smaller because smearing tends to make objects 
rounder, and therefore any quantity calculated from the 
cllipticities is smaller. The dispersion decreases because 
P 7 is much noisier than P s h- In the last two rows of the 



first block of Table |^ we also give the results for the cos- 
mic shear estimate if we apply neither PSF correction nor 
weighting. This demonstrates that the detection depends 
neither on the specific weighting scheme nor on the neces- 
sary smearing correction and can be seen in uncorrected 
data. However, in order to make an unbiased estimate of 
the shear dispersion, the PSF and smearing corrections 
must be taken into account and our particular choice of 
the weighting scheme minimizes the error on this estimate. 

In the second block of Table || we show the cosmic shear 
estimator if we apply the PSF corrections from the star 
fields individually. Note that the number of galaxies per 
field (and therefore the number of fields with N > N m i n ) 
also changes when correcting with different star fields, 
since the smearing correction (P s h/P sm )* is different from 
star field to star field (as noted in Sect. p~3| ). Therefore, 
for some galaxies P 7 is above or below the cut of 0.2 we 
introduced to reject objects with unphysically large el- 
lipticities. Comparing the various results for the cosmic 
shear measurement, we find (7 2 ) = (16.0 ± 1.2(PSF) ± 
11.8(stat)) x 10~ 4 for all fields, which shows that the dif- 
ference between the star field corrections is much smaller 
than the statistical error. 

In Fig. |l4| we plot the mean cosmic shear estimators of 
the galaxy fields when correcting with different star fields 
vs. number of galaxies per field. As noted above, due to the 
cut in P 7 the number of galaxies per field can differ. The 
very large vertical error bars are due to different numbers 
of galaxies per field (including galaxies with high elliptici- 
ties or not). However, if we remove this effect by applying 
the cut on the mean P 7 over all star fields rather than 
the P 7 obtained from correcting with an individual star 
field, we still find variations of the cosmic shear estimator 
per galaxy field up to A7 2 = 0.004, which is correlated 
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Table 3. Results for the cosmic shear for different PSF corrections. The first block shows the results if we do not 
correct for PSF effects: the first Col. indicates if we use raw ellipticities (e raw ) or anisotropy corrected ellipticities 
(e am ) and if we apply the smearing correction (P 7 ) or not (P s h)- The first row gives the fully corrected result [see 
Table ^] for reference. The next block shows the results when we apply PSF corrections from the individual star fields. 
The results quoted were obtained weighting individual galaxies [see Eq. (g)] with w = %n, requiring P 7 > 0.2 (or 
P s h > 0.2), and weighting the galaxy fields by Wi = N unless otherwise indicated. 
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with the mean shear per field. This can be understood by 
noting that 



1 



7 T , 



N n (N n - 1) 



(7n 
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(22) 



then A7 2 = 2 (j n ) A (7,1). The spread in the cosmic shear 
values obtained from different fields decreases for larger 
N, which shows that fields with more galaxies are less 
affected by noise. Note that the cosmic shear estimator 
can also be negative. 

5.4. Simulations 

In order to test the analysis pipeline we produced 140 
uncorrelated galaxy fields and 20 star fields (1' x 1'). 

Catalogues with intrinsic properties for the galax- 
ies were produced with the Stuff [] program (see 

1 Part of the Terapix software suite available at 
http : //terapix . iap . f r/sof t/ 



Erben et al. 2001). The galaxies were taken to be at the 
same redshift (z — 1). The effect of cosmic shear was simu- 
lated by doing ray tracing through a number of lens planes 
which were produ ced from N -body simulations in an open 
CDM cosmology ( Jain et al. 2000). Maps of the shear field 
produced from these simulations were sampled with ran- 
domly placed fields like the STIS fields. The cosmic shear 

signal on arcminute scales in these fields is J (t~ 2 ) = 2.6% 
( [Jain et al.| 2000). 

From the sheared galax y catalogues, images were pro- 



duced with Skymakern (see Erben et al. 2001), where the 



parameters were tuned to match as closely as possible 
the physical characteristics of the STIS CCD for exposure 
times of 400 s, so noise properties of individual simulated 
images are very close to the ones observed in real data. 
Star fields were created directly using Skymaker in the 
stellar field mode. We focused on obtaining stars of similar 
half-light radii as the observed ones, by setting the track- 



2 Part of the Terapix software suite available at 
http : //terapix . iap . f r/sof t/ 
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Fig. 14. For each of the 121 galaxy fields the cosmic shear 
was measured from ellipticities corrected with 21 individ- 
ual star fields. We show the mean of the cosmic shear 
estimator over the star fields as boxes vs. the number of 
galaxies in the field. The error bars show the la dispersion 
of individual values from the mean. Note that the number 
of galaxies per field can change because of the cut in P 7 . 



ing error type parameter in Skymaker to JITTER and its 
size to O'.'Ol. We simulated a round PSF and did not try 
to include the effects of optical distortion of HST, since 
the simulations were just done to test the analysis pipeline 
and its capability to retrieve the input cosmic shear signal. 
To simulate a coadded exposure of 1600 s, we first created 
4 individual 400 s images from the same original sheared 
galaxy catalogue, but with positions randomly dithered in 
order to process them in the same way as the real data. 
Dithers could vary between and 30 pixels. The pipeline 
was able to recover the dithering within one tenth of a 
pixel as stated in Paper I . The images were then drizzled 
and median averaged in the same way as the archive data. 

The analysis of the simulated fields was carried out fol- 
lowing the same procedure for catalogue production and 
cosmic shear estimation as described in Sect. ||and 5.1. In 
particular, we used the same parameters for simulated and 

real data. We find a cosmic shear signal of yj (^f 2 ) = 2.3% 
with 1.6a significance ((j 2 ) = (5.1 ± 3.2) x 10~ 4 ), which 



is slightly smaller than the input value of 2.6%. Since we 
under- rather than over-estimate the input signal we con- 
clude that our procedure does not introduce a spurious 
cosmic shear signal. 

For a quantitative analysis, we have to generate more 
realistic simulations which represent the properties of our 
data better with respect to different exposure times and 



therefore different redshift distributions, in order to better 
understand their contribution to the cosmic shear signal. 

6. Discussion 

We report on our initial measurement of cosmic shear 
on scales below one arcminute with STIS parallel archive 
data. Since any PSF anisotropy can mimic shear, we in- 
vestigated the PSF anisotropy of the STIS CCD in detail. 
Although the STIS PSF is not stable in time, we show that 
the anisotropy on the STIS field is sufficiently small (on 
the order of 1%) to carry out the cosmic shear analysis. It 
changes the mean shear on galaxy fields by much less than 
1% and does not affect our estimate of the cosmic shear 
rms by more than 10%. The smearing introduced by the 
PSF due to the small galaxy sizes at faint magnitudes is 
much more important than the PSF anisotropy and pro- 
duces the biggest dispersion in our results. This is a rather 
fundamental problem in that the galaxies are small and 
only NGST will have a higher spatial resolution. 

To test our procedure for catalogue production and 
cosmic shear estimation, we carried out simulations in- 
cluding the drizzling procedure. The cosmic shear result 
we obtained from the simulations is only slightly smaller 
than the input value, which shows that our procedure does 
not introduce a spurious cosmic shear signal and that our 
careful galaxy selection leads to a conservative estimation 
of the cosmic shear. 



As discussed in Sect. 5.2 weighting of individual galax- 
ies (according to the measurement error) or weighting of 
galaxy fields (to account for Poisson noise) changes the 
result substantially. If we restrict our analysis to fields 
with a higher number density of galaxies, we find that 
the signal increases on average, although not significantly. 
This would agree with a cosmological interpretation of the 
signal: fields with a higher number density of objects typi- 
cally have a larger exposure time, therefore they typically 
probe higher redshifts and one expects a higher cosmic 
shear signal. For a quantitative interpretation of the cos- 
mic shear measurement an optimal weighting scheme still 
has to be found, which will include redshift information 
on the STIS fields. 

Although we obtained a detection of the cosmic shear 
the error bars are still large. As can be seen in eqs. (12) 
and (|l|) the errors in the data depend on both the number 
of galaxies per field and on the number of fields. It is 
therefore important to get more fields with higher number 
densities of objects. The parallel observations with STIS 
are currently being continued in the frame of a GO cycle 
9 parallel proposal (8562+9248, PI: P. Schneider). 

In Fig. ^| we compare our result of 3.87t2;§4% for 
the cosmic shear to the ones obtained by other groups on 
larger scales and to the theoretically expected values when 
using different cosmological models with a mean redshift 
of the sources of (z s ) ss 1.2, which is appropriate for the 
groundbased measurements. With the STIS data we are 
probably probing at higher mean redshifts, but since our 
observations with STIS were taken in the CLEAR filter 
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Fig. 15. Comparison of our cosmic shear result if we use 
all galaxy fields (weighting individual galaxies by u>nn and 
galaxy fields with N) with measurements at larger angu- 
lar scales from other groups and with model predictions. 
9 is the radius in arcmin at which the results were ob- 
tained; the cosmic shear results which were measured on 
a square field (STIS, BRE, KWL) are plotted at the ra- 
dius of a circular field with the same area, the result from 
WTK is plotted at half the separation angle of the galax- 
ies. The lines show the theoretical predictions if one uses 
different cosmological models, which are characterized by 
fi m , 17a, h, r, and as- The redshift distribution is taken 
from Brainerd et al. (1996), with a mean source redshift 
of 



= 1.2. 



mode, we do not have much information about the redshift 
range of our galaxies. Moreover, our fields have a large 
spread in exposure times and therefore we effectively aver- 
age over different cosmic shear values. Multicolour (optical 
and near-IR) observations from the ground are presently 
being carried out to determine the redshift distribution 
of the galaxies in the STIS fields, using photometric red- 
shifts on the actual observed STIS fields and their vicinity. 
With these data, we will be able to interpret the cosmic 
shear measurement with respect to other scales. This will 
be done in a forthcoming paper. 
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